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SOLUTION  OF  A  MINIMUM  PROBLEM  IN  AIRCRAFT  WING  THEORY 


Karl  Nickel  in  Tubingen 


ABSTRACT 

The  third  fundamental  problem  of  Prandtl' s  theory 
of  the  supporting  line  is  extended  to  the  case  of  a 
finite  number  of  linear  supplementary  conditions,  and 
solved.  Three  examples  demonstrate  cases  which  occur 
in  practice,  with  such  supplementary  conditions. 

1 .  Formulation  of  the  problem 


In  wing  theory  I  of  the  "Four  Essays"1  \  L.  Prandtl  states  the 
following  minimum  problem  as  the  "Third  Fundamental  Problem"  (op. 
cit . ,  p .  28) : 

"Given  is  the  total  lift  and  the  wing  span,  also  p  and  V;  to  be 
found  is  that  distribution  of  the  lift  over  the  wing  span  for  which 
the  drab  becomes  a  minimum". 


With  the  aid  of  the  given  formulas  (op.  cit.,  p.  27)>  the 
mathematical  formulation  of  this  problem  reads  as  follows  in  a  coor¬ 
dinate  system  in  which  the  wing  extends  from  -1  to  +1: 


Let 


1 

4*_jT  iy  * — ( 


(n 


_  und  Aerodynamik 

Prandtl,  L.,  and  Betz.,  A.  Vier  Abhandlungen  zur  Hydrodynamlk^ 
Neudruck  aus  den  Verhandlungen  des  HIT  Internatlonalen  Mathem- 
atlker-Kongresses  zu  Heidelberg  und  aus  den  Nachrlchten  der 
Gesellschaft  der  Wlssenschaften  zu  Goettingen^  (""Four- Essays 
on  Hydrodynamics  and  Aerodynamics.  Reprint  from  the  Proceedings 
of  the  Third  International  Congress  of  Mathematicians  at  Heidelberg 
and  from  the  Communications  of  the  Society  of  Sciences  at 
Goettingen").  Goettingen,  1927. 

P  \ 

'  The  sign  is  to  be  understood  as  the  Cauchy  principal  value  of 
the  integral. 
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and  _ i) _ _ 0  [where  r(x)  is  the  local  lift  density  and 

w(x)  the  induced  downwash  velocity  at  position  x  on  the  wing].  Let 
the  function  T(x)  on  the  interval  <-t  +1>  be  determined  such  that 

}  r(x)\c(x)dx  C2> 

~1 

becomes  a  minimum  under  the  supplementary  condition 

j'r(x)dx=A  (3) 

(A  =  the  total  lift,  arbitrarily  prescribed). 

■3  ) 

This  formulation  of  the  problem  was  extended  by  M.  Munk  to 
the  case  of  arbitrarily  distributed  and  directed  lift,  and  solved 
in  general. 

One  can  also  generalize  the  above  minimum  problem  in  another 
way,  by  prescribing,  instead  of  the  supplementary  condition  (3), 
other — and  if  necessary  several — supplementary  conditions.  Three 
examples  will  illustrate  this  more  exactly: 

a)  A  flying  wing  is  to  be  made  into  a  (shallow)  curve  by  a 
(small)  aileron  deflection.  What  form  must  be  taken  by  the  lift 
distribution  which  is  added  to  the  symmetric  distribution  by  the 
aileron  deflection  in  order  that  the  increase  in  the  induced  drag 
stays  as  small  as  possible?  If  one  knows  this  distribution,  one  can 
approximate  it  with  a  suitable  aieleron  shape,  in  order  to  prevent 
superfluous  losses  when  deflecting  the  ailerons.  The  question 
here  therefore  reads:  The  lift  distribution  is  sought  which  makes 
the  induced  drag  (2)  a  minimum,  when  the  roll  torque 

Jr(x)xix  (V) 

is  to  have  a  prescribed  value.  The  use  of  equation  (1)  is  in  this 
case  exactly  correct  only  in  the  first  instant  of  aileron  deflection, 
so  long  as  there  is  no  yawing  or  rolling  motion,  because  only  then 
the  departing  vortex  lines  are  rectilinear.  At  any  rate,  for  small 

T5 - 

Munk,  M.  Isoperlmetrlsche  Aufgaben  aus  der  Theorle  des  Fluges . 
Inaugural-Pis  sert  at  lorTI  C"Isoperimetric  Problems  from  the 
Theory  of  Flight.  Inaugural  Dissertation").  Goettingen,  1919. 
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yaw  and  roll  velocities  one  can  still  use  the  solution  of  the  for¬ 
mulated  problem  as  an  approximation. 

A  second  problem,  which  leads  to  the  same  (mathematical)  prob¬ 
lem,  and  in  which  these  difficulties  do  not  arise,  was  kindly  commu¬ 
nicated  to  me  by  Dr.  Prandtl.  This  is  the  question  of  the  lift 
distribution  with  minimum  drag,  for  a  wing  with  an  eccentrically 
applied  load,  which  creates  a  roll  torque  (A). 

b)  The  solution  of  the  third  fundamental  problem  stated  by 
Prandtl  reads  (ref.  1,  page  32)  in  the  formulation  given  above: 

that  is,  the  lift  is  distributed  over  the  wing  span  in  the  form  of 
a  half  ellipse  (see  the  solid  line  in  Figure  1).  If  one  now  modi¬ 
fies  T(x)  somewhat,  in  the  manner  shown  by  the  dashed  line,  letting 
A  in  (3)  keep  its  value,  the  induced  drag  (in  the  vicinity  of  its 
minimum)  will  change  only  slightly.  The  spar-bending  torque  at  the 
wing  root 

j  r(x)xdx, 

w'hich  can  be  replaced  by 

Jr(x)\x\dx  (y) 

for  symmetric  lift  distributions  (apart  from  a  factor  of  1/2),  will 
however  become  smaller. 

For  a  free-flying  wing,  this  means  the  following:  in  conse¬ 
quence  of  the  reduced  stress,  the  wing  can  be  built  lighter  in 
weight.  Because,  however,  in  constant  horizontal  flight  the  lift 
and  the  weight  are  equal,  this  means  that  the  total  lift  A  in  (3) 
also  becomes  smaller. 

If  F(x)  is  replaced  by  X*T(x),  the  lift  (3)  and  spar-bending 
torque  (5)  change  in  proportion  to  X,  while  the  induced  drag  is  pro¬ 
portional  to  X2.  One  would  expect,  therefore,  somewhat  more  favor¬ 
able  drag  conditions  with  a  lift  distribution  which  has  somewhat 
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Figure  1  Figure  2 

smaller  values  at  the  wing  tips,  than  with  the  elliptical  lift  dis¬ 
tribution.  At  any  rate,  it  must  be  noted  that  the  dependence  of 
the  flying  weight  on  the  spar  bending  torque  is  very  small  for 
the  usual  wing  construction  methods,  so  that  the  effect  described 
above  is  small. 

This  consideration  leads  to  the  following  formulation  of  the 
problem:  The  induced  drag  (2)  is  to  be  made  a  minimum  under  the 

two  supplementary  conditions,  that  the  lift  (3)  and  the  spar  bending 
moment  (5)  possess  prescribed  values. 

c)  If  one  considers  a  wing  which  is  curved  in  the  direction 

of  flight  (see  Figure  2),  the  longitudinal  torque  would  also  have 

to  be  held  constant,  in  addition  to  the  total  lift,  in  solving  the 

in 

third  fundamental  problem  of  wing  theory  for  this  case.  One  can 
now  require  that  the  constancy  of  the  longitudinal  torque  should 
result  exclusively  from  lifting  forces  on  the  wing  (longitudinal 
stabilization  through  wing  warping,  such  as  in  bird-flight).  Let 
h(x)  be  the  distance  of  the  central  line  of  pressure  (dashed  in 
Figure  2)  from  the  x-axis.  Then  this  requirement  means  that 

//X*)  *(•)'*«*  (<*> 

—  I 

51 - 

It  is  true  (see  reference  1,  page  25)  that  the  idea  of  the  sup¬ 
porting  line  can  no  longer  be  used  in  this  case;  however,  as  M. 
Munk  has  shown  (reference  3»  page  21),  in  considering  the  drag  it 
suffices  to  consider  only  the  two-dimensional  problem  of  a  wing 
which  is  not  curved  in  the  direction  of  flight.  T(x)  here  is 
thus  the  projection  of  the  lift  onto  a  plane  normal  to  the 
direction  of  flight. 
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(longitudinal  torque  M  prescribed)  and  the  problem  above  reads: 

The  lift  distribution  r(x)  is  sought,  for  which  (2)  becomes  a  min¬ 
imum  under  the  two  supplementary  conditions  (3)  and  (6).  It  might 
be  remarked  here  that  in  the  general  case  of  an  arbitrarily  shaped 
wing  the  central  line  of  pressure  will  not  be  the  line  of  the  for¬ 
ward  neutral  point  in  the  two-dimensional  problem  (T/4  line).  The 
task  of  determining  this  line  for  an  arbitrarily  prescribed  wing 
poses  great  difficulties  which,  to  my  knowledge,  have  not  yet  been 
solved.  As  an  approximation,  one  can  nevertheless  replace  the  cen¬ 
tral  line  of  pressure  by  the  T/-U  line,  so  long  as  the  wing  does  not 
deviate  too  much  from  a  straight  wing. 

On  the  basis  of  these  examples,  it  is  now  natural  to  extend 
the  minimum  problem  formulated  above  in  such  a  way  that  finitely 
many  supplementary  conditions,  linear  in  T(x),  are  prescribed.  One 
seeks  therefore  to  determine  T(x)  in  the  interval  < — +1>  such  that 

jr(z)v(z)dz  (2s) 

—i 

becomes  a  minimum,  under  the  N  supplementary  conditions 

jf/x*)  ».(*)<*«  4.  (*=t . N) 

[JX— l)=r(+l)  =  0;  *,(«)  and  An  arbitrarily  prescribed,  except  that  for 
h  (x)  =  0  naturally  only  A  =  0  is  allowed].  For  N  =  l;  *,(*)  =  l,  4,  =  A 
the  problem  formulated  by  Prandtl  results  as  a  special  case;  like¬ 
wise  the  examples  a),  b)  and  c)  can  be  obtained  through  analogous 
specialization . 

2 .  Transformation  of  the  problem 

For  the  mathematical  treatment  of  the  problem,  let  the  follow¬ 
ing  transformations  be  undertaken: 

One  sets 


•/<«). 


far  CM  f, 
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and  uses  the  notation  *»(x) /l  —  x*  =  *,(j)  (*  =  1.  . . N)  •  Then-5'  from 
(1),  (2)  (neglecting  a  factor  of  1/4)  and  (7): 


JW  *7  it  cos  t  — —  cos  $ 

0 


~^JZ(i)f(i)ds 


(»■). 

(2a), 


(«=1 . N)  ....(7a). 


For  known  f(s),  (la)  is  a  Fredholm  integral  equation  of  the  first 
kind  for  dZ(t)/dt.  As  is  well  known,  its  solution^  is 


ijLLi 

it  x 


sin  tit 
cos  9  —  cos  t 


from  which  follows  by  integration 

* 

Z{»)  I  AQ  log 


.  *+* 

sin  — sr~ 


sin 


«  — I, 


it 


if  one  writes  for  brevity 


.  «  +  * 
1,  ,W-2 
*  °®  .  !•— I, 

“nr 


5  (*.«). 


and  inserts  (8)  in  (2a)  and  (7a),  one  obtains 

///(*)/(«)  8(*,t)d,dt 


and 


j  j A*)  t)dsdt  =  AH 


(V 


(2  ><o 

(j£) 


Since  the  functions  Z(s)  and  f(s)  are  unambiguously  related  to 


5 ) 

The  conditions  under  which  the  integrals  and  infinite  series  which 
occur  exist,  and  under  which  the  transformations  used  (changing 
the  order  of  integration,  exchanging  integrations  and  summations) 
are  allowed,  are  examined  in  a  comprehensive  work,  which  is  to 
appear  shortly  under  the  title  Solution  of  a  Particular  Minimum 
Problem  in  the  Mathematische  Zeitschrift.  In  it  a  second  class 
of  supplementary  conditions  will  be  added  to  those  already  consi¬ 
dered.  Note  added  in  proof:  This  has  meanwhile  appeared  in  the 
Mathematische  Zeitschrift,  vol.  53  (1950),  pp.  21-52. 

^See  for  example  Schroeder,  K.  An  Integral  Equation  of  the  First 
Kind  in  Wing  Theory.  Sitzungsberichte  der  Preussischen  Akademie 
der  Wissenschaften  XXX  (1938). 
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each  other  through  (la)  and  (8),  one  can  also  consider  f(s)  as  the 
function  to  be  sought.  With  this  the  problem  to  be  solved  reads: 
f(s)  is  to  be  determined  in  the  interval  <0,*>  ,  such  that  (2b) 
becomes  a  minimum  under  the  supplementary  conditions  (7b). 


For  further  simplification  of  notation,  let  the  following 
further  abbreviation  be  introduced: 

(?) 


«  0 


With  the  representation 


7) 


B{t,  t)  =  —  tit  sin  «t 


valid  for  0  < », t  t  +  i  ,  there  results  from  (9): 


(/.*)  =  /  //(«)j(‘)S(M)i«  dt 

*  *  -  »  \f*  \ 

=  j-  -i  |  j /(*) sin  ns  dsj M  g(t)  sin  ntitj, 


(?o) 


with  the  Fourier  sine  coefficients 

2  f  It 

•n~-J  /(•)*>»•*•,  b.  =  -J  ,(t)»inntit  (n=l,  2,...). 

•  • 

According  to  (9)  and  (10),  and  because  S(s,t)  =  S(t,s),  (f,g)  has 
the  following  properties:  (a  =  a  real  number) 


*)  («/,  g)  =  « (J,  g)  (a  »=  reelle  Z»hJ) 

b)  ,,  (/•?)  =  >/) 

c)  </+?.*)= a. *) +(?.*) 

<*)  (/,/)>  0  fflr  /(*)  -:|  0 


r>/> 


One  therefore  has  the  task  of  finding  a  function  f(s)  in  the 
interval  <0,»>  such  that  (f,f)  becomes  a  minimum  under  the  N 


7^See  for  examp.^.  H:  si,  G.  Integralglelchungen  ("Integral 
Equations"),  u-rlin:  J.  Springer,  1937  ;  or  Jaeckel,  K. 
Determination  of  a  Series  Representation  for  the  Kernel  in  r 
in  Elliptical  Coordinates.  Z.  angew.  Math.  Mech.,  vol.  30 
(1950),  p.  186  formula  (16) . 
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(/»/)  =  Min. 


supplementary  conditions 

(/.*.)=  **  (»=1 . N)  (/£> 

Thereby  (f,g)  is  defined  by  (9)  and  has  the  properties  (11). 


3 .  Solution  of  the  problem 

If  the  functions  h  (s)  are  linearly  dependent,  then  under  the 
supplementary  conditions  (12)  either  they  are  supernumerary  or  they 
are  contradictory  and,  therefore,  cannot  be  satisfied  at  all. 
Therefore,  without  loss  of  generality,  one  may  assume  the  functions 
hn(s)  to  be  linearly  independent. 

Assumption:  Let  the  functions  h^(s)  be  linearly  independent, 

and  let  the  real  numbers  A^(n  =  1,...,  N)  be  arbitrarily  chosen. 

Proposition:  Then  there  is  exactly  one  solution  f!(s)  to  the 

above  minimum  problem  and  this  solution  has  the  form 

/'(')= 2  «.*.(«) 

«  ! 

with  unambiguously  determined  real  numbers  a  . 


Proof:  According  to  the  orthonormalization  procedure  of  E. 

Schmidt,  one  can  find  functions 

Hn($)  —  2  e mnkm(t)  (c.,  +  0;  »  =  1 . N) 

mml 

such  that 


II  for 

(Hm,  H.)  =  «***  =  |0  for 


ni*n 
*1  +  " 


Inverting  these  equations  gives  the  representation 


KM -JZKnB-M  (Kn  +  0;  n -  1.  . . N). 


One  now  sets 


/'(•)*  i 

••  i 

with  real  coefficients  c^  whic_.  are  as  yet  undetermined  and  with 
the  calculation  rules  (11)  obtains 
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(/'.  K)  =  2  ‘n(Hm,  h,)  =  2  cm  2  imr(Hm,  Hm) 

*— l  »-l  ■*— l 

'2*n  2imr**m*=  f0r  f=l . N. 


■  -1 


r 

If  now  2  cn&nf  *=  -4,  (f  =  1,  . . tf)  ,  the  are  thereby  unambiguously 

determined  in  a  recursive  fashion.  With  the  real  numbers  c  thus 

n 

chosen,  ff  then  satisfies  the  supplementary  conditions  (12),  and 
further  receives  the  desired  form 

•/'(#)  =  i  CW Hm (4)  =  i  c.i  A* M«)  =  i  M*)  i  e-  =£  °-  *»  (*>• 

a-1  a— l  K-l  »-l  •— •  •"* 


It  still  remains  to  show  that  f*  also  minimizes  (f,f).  An 
arbitrary  function  f(s)  can  always  be  written  in  the  form 

f(s)  =  f'(s)  +  k(s)  (namely  with  k(s)  =  f(s)  -  f'(s)). 

Thus 

(/A)  =  (/'  +  *■*.)  =  (/A) +(*.*.)==  4 .  +  (*A)  (n=l . N). 

f(s),  therefore,  satisfies  the  supplementary  conditions  (12)  exactly 
when  (k,hn)  =  0  for  n  =  1,...,N. 

For  all  the  functions  f(s)  which  satisfy  the  supplementary  con¬ 
ditions  (12),  the  following  holds,  according  to  the  rules  (11): 

(/./)  =(/'+*•  /'  +  *)  =  (/'./') +  2  (f. k)  +  (*.  k) 

=  (/'./')  + 2 1  *)+(*,*) 

whereby  the  equality  sign  comes  into  play  only  for  k(s)  =  0.  The 
unambiguity  of  the  coefficients  an  follows  with  (lid)  for  the  linear 
independence  of  the  functions  h  (s).  Thus,  the  proposition  is  proven. 

If  one  reverses  the  transformations  introduced  above,  one  finds 
the  following  result: 

4.  Result 


Let  the  functions  hn(x)  and  the  real  numbers  A^  (n  =  1,...,N) 
be  arbitrarily  chosen.  The  function  r(x)  in  — is  to  be 
determined  such  that  and  such  that  by  setting 
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the  expression 


jV(*)«(*)i*  (z  ^ 

becomes  a  minimum  under  the  supplementary  conditions 

Jr(*>  *»(*)<*»  (««i . N)  (7) 


If  a  function  T(x)  exists  at  all,  which  satisfies  the  N  supple¬ 
mentary  conditions  (7),  then  this  minimum  problem  has  exactly  one 
solution,  which  has  the  form 


i 


—  *»  +  1  —  ^)U  —  I*)! 


with  real  coefficients  an .  If  the  functions  h  (x)  are  linearly 
independent,  then  the  minimum  problem  always  has  a  solution  and  the 
coefficients  an  are  unambiguously  determined. 


5 .  Application 


In  the  following  table  are  the  associated  functions  Mx) 


Func t ion 


/  *(y)  iog  y*)[  iy 

-i  i  *— »  4  y‘ 


fw..x,.g|i-*v  +  V(i -**)(!-•/) 


**  g-(2**+  1)>1  —  ** 

I  1*>»  2l,0|!l+t'!p? 

I— 1*<0  R  1*1 

1*1  **,0«L+|I| |— 

,*l*l  !(**log1+^|  ■**  +*Vl  — **) 

**ixi  1  (e  log 1  *  £ + (6 »» + 1)  yr=?) 

(a» -x)iog|1~a,g+ -1- — 1 + }fvzr* arc cos *• 

With  these,  one  finds  the  solutions  for  the  first  two  cf 
examples  stated  previously: 

x)  jl\x)xix^  R  r(x)  =  ^ 


jr/x*)*i*=»  r(x)=  h  *ii-x» 

|  jr(x)dx=A  r(x) = s  — h) 

b)  I  »~‘  *_  .  44\  .... 


i+ym? 
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In  Figure  3  is 
4  shc.is  some  curves 
for  A  =  constant,  H 
4A  the  half-ellipse 


displayed  this  solution  for  example  a).  Figure 
from  the  family  of  solutions  for  example  b) 
variable  (among  these  is  contained  for  3tH  = 
as  a  special  case). 


Figure  4 


Submitted  on  January  28,  1950. 
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